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Abstract 

P 

■ Variational calculations of the deuteron and the triton illustrate that simple wave 

function ansatze become more effective after evolving the nucleon-nucleon poten- 
tial to lower momentum ("Viowfc")- This is consistent with many-body wave func- 
tions becoming much less correlated at lower cutoffs, as shown by two-particle wave 
functions and pair-distribution functions in nuclear matter. These results motivate 
| a program to explore variational many-body calculations of binding energies and 

' other low-energy nuclear properties using low-momentum potentials. 
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1 Introduction 



Variational many-body calculations rely on being able to devise and optimize 
sufficiently rich wave function ansatze. The strong short-range repulsion and 
tensor forces of conventional nucleon-nucleon potentials that fit phase-shift 
data necessitate complicated, highly correlated trial wave functions, which 
limits the effectiveness of the variational approach. Even large-scale calcula- 
tions of light nuclei using variational Monte Carlo (VMC) achieve accuracy 
only at the five percent level [1] and Green's function Monte Carlo (GFMC) 
is needed to evolve the wave functions [2] . Recent results with low- momentum 
potentials, however, suggest that more effective variational calculations should 
be possible [3]. 

The nonperturbative nature of conventional inter-nucleon interactions can be 
radically modified by using the renormalization group to lower the momentum 
cutoff of a two-nucleon potential [4]. For low-momentum interactions with 
cutoffs around 2fm _1 , the softened potential combined with Pauli blocking 
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leads to corrections in nuclear matter in the particle-particle channel that 
are well converged at second order in the potential. Calculations of nuclear 
matter using the low-momentum two-nucleon force V\ ow k with a corresponding 
leading-order three-nucleon (3N) force from chiral effective field theory (EFT) 
exhibit nuclear binding in the Hartree-Fock approximation, and become much 
less cutoff dependent upon adding the dominant second-order contributions. 
At the lower cutoffs, the iterated tensor interaction in the two-body sector does 
not play a major role in nuclear saturation, in contrast to conventional wisdom. 
Thus, the relative importance of contributions to observables from the tensor 
force or from three-body forces is strongly scale or resolution dependent. 

Similarly, the correlations in many-body wave functions will exhibit signif- 
icant resolution dependence. The dominance of Hartree-Fock and the onset 
of perturbative convergence in the particle-particle channel at smaller cutoffs 
implies that the corresponding wave functions are much less correlated than 
those associated with conventional potentials. This has the practical conse- 
quence that variational calculations should be effective with much simpler 
ansatze. In this letter, we illustrate this simplicity through variational cal- 
culations of the deuteron and triton, and by examining the pair-distribution 
and two-particle wave functions in nuclear matter at empirical saturation den- 
sity. We are not trying at this stage to optimize the variational approach for 
low-momentum potentials. Rather our goal is to motivate a program to re- 
examine the application of variational methods to binding energies and other 
low-energy properties of nuclei using low-momentum interactions. 

The construction of V\ ow k is well documented in Refs. [5,6], where it is shown 
that either renormalization-group equations or Lee-Suzuki transformations 
can be used. Here we employ the latter, which provide a convenient formal- 
ism to evolve consistent operators beyond the Hamiltonian. In the notation 
of Ref. [7,8], the evolution of an operator O from a momentum cutoff A to 
A < A is given by 

0(A) = * , (P + ^)d(A )(P + u) ] , (1) 

where the operator oj = QujP parameterizes the Lee-Suzuki transformation, 
the projection operator P projects onto relative momenta k < A, and Q 
projects onto A < k < Aq. In the case of the evolved Hamiltonian it is con- 
venient to define V[ ow £:(A) = H(A) — T, where T is the "bare" kinetic energy 
operator. By construction, two-body bound-state properties and phaseshifts 
are preserved for external relative momenta up to the cutoff. Three- and many- 
body observables require the consistent addition of a three-body force to re- 
move cutoff dependence [6]. 

We will show results starting from the Argonne Vi$ potential [9], since it is 
used in almost all modern VMC and GFMC calculations. However, for cutoffs 
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of 2fm _1 or less, all bare potentials that reproduce nucleon-nucleon phase 
shifts up to 350 MeV lab energy, including EFT potentials at N3LO, collapse 
to the same V\ ow k [5]. Therefore, the pattern of results for the full cutoff range 
shown here does not vary significantly with different initial potentials. 

Since Vj ow fc is energy independent, variational calculations with V\ ow k proceed 
as described in ordinary quantum mechanics texts (e.g., without special nor- 
malizations as needed for energy-dependent potentials). That is, given a trial 
wave function Atrial) our variational estimate for the ground state energy at 
cutoff A is: 

jp ( \ \ (Atrial 1^+ Mow fc(A) Atrial) (e) , 
-C/varlAJ = — — r , (2 

(Atrial | Atrial) 

which we minimize with respect to the parameters in Atrial- Alternatively, we 
get a variational estimate by diagonalizing T + Vj ow fc(A.) in a truncated basis, 
where the trial wavefunction is a linear combination of the basis functions. 



2 Deuteron 



We start with a study of the deuteron binding energy. The philosophy is that 
for a weakly bound state one shouldn't need to work hard, so a simple, generic 
ansatz should work well. We test this as a function of the cutoff. Our first 
ansatz is adapted from the discussion of wave functions in momentum space 
given long ago by Salpeter [10], which motivates the (unnormalized) 3 Si and 
3 £>! trial functions for the deuteron (following the conventions of Ref. [11]) 

^) = (Jfc 2 +7 2 )(A .2 + /i2) > Mk) = (Jfe2 ~ ^(p + ^2 » ( 3 ) 

where 7, fx, v, and a are variational parameters. Obviously we could extend 
this ansatz in many ways, but our point is to see how well we can do without 
having detailed knowledge about the wave function. The underlying physics 
implies that \x and v should be roughly the inverse range of the interaction 
and 7 should be close to (— M^E^I 2 , where M N is the mean neutron-proton 
mass and E d pa —2.2246 MeV is the deuteron binding energy. Moreover, the 
cutoff in Viowfc implies that the exact deuteron wave function does not contain 
high-momentum components. Therefore, the two-body trial wave functions 
considered contain the same cutoff on the relative momentum. 

The best variational energy for Eq. (3) as a function of the cutoff is shown as 
the squares in Fig. 1. These estimates are not even bound for cutoffs above A pa 
5 fm _1 (which includes the bare Argonne v±s potential) but rapidly improve as 
the cutoff is lowered further, reaching a minimum deviation of less than 3 keV 
around A pa 1.5 fin -1 . We emphasize that the low-momentum potential really 
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Fig. 1. Deviation from Ed of the best variational energy as a function of cutoff A 
for the wave function ansatze of Eqs. (3) and (4) with different numbers of terms. 

does preserve the two-body observables; directly solving the momentum space 
Schrodinger equation with 40 gauss points yields the same deuteron binding 
energy as the bare potential to a fraction of an electron volt for all the cutoffs 
shown. But variational calculations with the lower cutoffs come much closer 
to this energy. 

We can also adapt the form used for a high-accuracy representation of the 
deuteron wave function in Ref. [11] (and elsewhere) to see if the same pattern 
holds. Consider 

n ft n 7~) 

~{ k l + rrij jr[ Ar + rrij 

where the rrij are fixed in a geometric progression: 

m . = {-M N E d f /2 + (j - l)m , with m = 0.9 far 1 , (5) 

and we treat the Cj and Dj as variational parameters for a given value of 
n. (The very accurate parameterization of the deuteron wave function for the 
Bonn potential in Ref. [11] has n — 11 with some constraints on the Cj's 
and Dj's.) Since the variational coefficients appear linearly, we can simply 
diagonalize the Hamiltonian in the truncated basis of Eq. (4) to find the best 
estimate of the deuteron energy. For n = 2, the best estimate is about 5keV 
above Ed at A « 1.4 fm _1 and is unbound above 6 fm _1 . By enlarging the basis 
(n = 3 and n = 4 are shown in Fig. 1) we find better estimates for all A, but 
the steep improvement remains with the minimum shifting gradually higher. 
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The success of simple ansatze for the deuteron at low cutoffs can be under- 
stood by looking at the corresponding wave functions. In Fig. 2, we show the 
exact deuteron wave functions in momentum space for a variety of cutoffs. The 
immediate source of the problem for a good energy estimate at higher cutoffs 
is the node in momentum space, which reflects the short-range correlations 
that are evident as a suppression of the coordinate-space wave function for 
r < 1 fm _1 (see Fig. 3) [12]. This short-range/high- momentum behavior is in- 
creasingly resolved at higher cutoffs, which entails finer and finer cancellations 
in the variational integrals. 

As suggested by the plots of the 3 Si component ipo(k) in Fig. 2 and the anal- 
ogous behavior of the 3 Di component tp2(k), matrix elements of the operator 
a k a k (which are proportional to |?/>o(^)| 2 + |"02(^)| 2 ) change as A is lowered. 
This is neither surprising nor worrisome, since the momentum distribution 
defined by the "bare" a\ i a k operator is not an observable (see Ref. [13] for a 
detailed discussion), and more generally because operators beyond the Hamil- 
tonian must also be evolved in order to give A independent expectation values 
[14] . However, as a test of our crude variational wave functions, it is instructive 
to start with this operator and a potential at a large cutoff A (e.g., 10 fur 1 ) 
and evolve both down to A using Eq. (1). By construction, matrix elements 
with the exact wave functions as in Fig. 2 are unchanged even for k > A; 
we have verified this explicitly. One might worry that an intricate interplay 
of operator and wavefunction is needed to preserve the matrix element and 
this might be lost with a variational wave function, but this is not the case 
with low momentum interactions. For example, matrix elements of the evolved 
operator with the simple variational wave function of Eq. (3) for A 1.5 fm -1 
(which gives the best energy estimate) underpredict the exact result by less 
than 4% for k < 0.1 fm -1 and then closely reproduce it for all k up to 5fm -1 . 
This issue will be examined in more generality in future work. 



3 Triton 



Moving from the deuteron to the triton is a significant step in complexity, but 
we can still test fairly simply whether the basic deuteron results carry over by 
using a truncated harmonic oscillator basis for a variational calculation. The 
antisymmetric three-nucleon basis is generated from the Jacobi coordinate 
oscillator states [8] 

\(rdsjt;Mc\j\)JT), (6) 

where (nlsjt) and {M C\j\) are the quantum numbers corresponding to the 
two relative Jacobi coordinates [e.g., k = — P2) and q = |(p 3 — \(pi + 
p 2 ))], and the basis is truncated according to the total number of oscillator 
quanta N = (2n+l+Z\f+C) < N max . Diagonalizing the intrinsic Hamiltonian 
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Fig. 2. Momentum-space 3 Si deuteron wave function for a range of cutoffs and for 
the bare Argonne vis potential. 
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Fig. 3. Coordinate-space 3 Si deuteron wave function for a range of cutoffs and for 
the bare Argonne v\$ potential. 
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Fig. 4. Relative deviation of the variational energy from the Faddeev result for the 
triton (Et) with two-body interactions only, as a function of cutoff A for several 
truncated oscillator basis sets. 



in the truncated basis and minimizing with respect to the oscillator length 
parameter b provides a variational estimate to the true ground-state energy. 



Results for the triton with two-body interactions only are shown in Fig. 4. 
Since the three-body contribution varies with the cutoff, the reference value 
for two-body alone varies as well. Therefore, the results are given at each A 
with respect to the exact Faddeev result [6] at that A, which we label E t . 

We see in the figure the same qualitative pattern as found for the deuteron. We 
can calibrate the efficacy of the harmonic oscillator basis at lower cutoffs by 
comparison to the calculation by Nunberg et al. [15], which used the Reid soft- 
core interaction (which actually has quite a repulsive core). Even using a basis 
with an additional nonlinear variational parameter compared to the current 
calculation, their predicted triton energy is not even negative until A^ max > 12 
and the largest calculation with A^ max = 28 yields E t = — 6.7MeV, which is 
extrapolated to E t = — 7.3±0.2 MeV for A^ max — > oo. That is, A^ max = 28 is still 
600 keV above the converged result. In contrast, for Vi ow fc with A = 1.8 fm -1 , 
the result for A^ max = 6 is already only 620 keV above the converged result, 
and this drops to 160 keV for A" max = 20. 1 



1 We anticipate even better variational convergence properties if one uses a smooth 
momentum cutoff in Vj ow fc- 
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4 Nuclear Matter 



We are optimistic that the general pattern we have seen in simple two- and 
three-nucleon variational calculations will continue for heavier systems (and 
with the inclusion of the corresponding low-momentum three-body force), and 
that the variational improvement upon lowering the cutoff is universal. Our 
optimism is based on calculations of nuclear matter with low-momentum po- 
tentials, which are discussed in Ref. [4]. 

Two convenient measures of correlations in nuclear matter are the in-medium 
pair wave functions and the pair-distribution function. Working within the 
Brueckner approximation, the relative in-medium wave function for a pair of 
nucleons with total momentum P and relative momentum k is 

l^)=|k) + -^7rG P M|k), (7) 

where Q is the Pauli blocking operator and G is the usual Brueckner G matrix, 
which sums the in-medium particle-particle ladder diagrams. In the current 
calculation, all self-energy effects are neglected for simplicity, which corre- 
sponds to setting u = k 2 + ^P 2 and using the pair kinetic energy operator for 
Hq. We stress that the overall pattern of our results are not changed by this 
assumption. At the same level of approximation, the pair-distribution function 
is given in terms of the coordinate-space pair wave functions [16], 

s(r)=X;(2T+l)(2S+l) £ \*Z(r;ST)\ 2 . (8) 
ST lf±k| 

Physically, g(r) corresponds to the correlation function for finding another 
nucleon a distance r from the first. 

In Fig. 5 we show the 3 Si wave function in the 3 S±- 3 Di coupled-channel for 
a pair of nucleons with P = and k = O.lfm -1 , and in Fig. 6 the pair- 
distribution function in nuclear matter at empirical saturation density k F = 
1.35 fm _1 . Unless otherwise noted, the curves include the dominant effects 
of the corresponding three-body force at each cutoff by converting the 3N 
vertex into a density-dependent two-body correction that is added to Vi ow fc, 
see Ref. [4] for details. 

The two-particle wave functions and pair-distribution functions in nuclear 
matter exhibit the same promising features we found in the simpler two- 
and three-nucleon systems. Namely, the strong short-ranged correlations are 
"blurred out" as the interactions are evolved to lower momenta. This is clear 
for the two-particle wave function, where the "wound" resulting from the 
short-range repulsion is more pronounced at the higher cutoffs. Similar re- 
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Fig. 5. Two-particle wave function for the 3 Si channel in symmetric nuclear matter 
at = 1.35 fm . 
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Fig. 6. Pair-distribution function in symmetric nuclear matter at kp = 1.35 fm - . 
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suits are found for the pair-distribution function, where g(r) at smaller cutoffs 
has little short distance structure and lies fairly close to the Fermi gas (i.e., 
Hartree-Fock wave functions) values, where the correlations arise solely from 
Fermi statistics. 

It is interesting to note that the correlations induced by the three-body force 
are significantly stronger at larger cutoffs. At A = 3.0 fm -1 , which is the 
largest cutoff at which the three-body force was fitted for Vi ow fc [6], the pair- 
distribution function is suppressed at short distances (relative to the Fermi 
gas values) when the three-body force is included. In contrast, there is a slight 
enhancement in g(r) at short distances when only the two-body V\ ow k is used. 

One finds that the changes in g(r) with and without the three-body force are 
much less severe at lower cutoffs. For example, at a cutoff of A = 1.9 fm -1 
one finds that g(r = 0) decreases by only 0.1 when the three-body force is 
included, while the corresponding change for A = 3.0 fm -1 is a factor of two 
larger. This is consistent with the results of Ref. [6], where it was shown that at 
smaller cutoffs the three-body force can be accurately treated in perturbation 
theory for the triton and alpha particle ground state energies. 

As we saw for the deuteron, these plots show that short-range correlations do 
not need to be explicit in the wave functions when calculating low-energy/low- 
momentum observables up to nuclear matter densities. The key point is that 
lowering the resolution makes calculations simpler, more efficient, and less 
mo del- dep endent . 



5 Conclusions 

In summary, the direct evidence from the deuteron and triton, coupled with the 
rapid convergence of the particle-particle channel observed in nuclear matter, 
imply that low-momentum potentials with A 2 fm~ will be much more ef- 
fective for few-body and many-body variational calculations than any conven- 
tional large-cutoff potential. Furthermore, even chiral EFT potentials, which 
are themselves low-momentum potentials compared to conventional potentials 
such as Bonn, Nijmegen, or Argonne, can be made more effective by running 
the cutoff lower. The general idea is to take the EFT cutoff as large as possible 
(i.e., in the vicinity of the breakdown scale of the chiral EFT, which presum- 
ably is 3fm _1 or higher), in order to minimize the truncation error [17]. The 
evolution to lower cutoffs induces the higher-order short-ranged operators that 
maintain the same truncation error in observables as at the higher cutoff. 

To take advantage of these observations, variational Monte Carlo (VMC) is 
attractive for its basic simplicity, the upper-bound property of the energy es- 
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timates, and the absence of the fermion sign problem. Variational calculations 
with conventional (large cutoff) potentials are performed in coordinate space, 
where the strong correlations are most naturally encoded in trial wave func- 
tions [1]. In fact, most calculations of this sort have used some variety of the 
Argonne potential, which was designed for this purpose and for subsequent 
GFMC calculations built on the variational results (local and operator based) 
[1,2]- 

Based on our results here, we anticipate more efficient variational results for 
low-momentum interactions, with the added advantage of being able to vary 
the cutoff as a tool to optimize and probe the quality of the solution. We have 
no restriction to coordinate space and, indeed, we will first try to develop 
Monte Carlo calculations for light nuclei directly in momentum space. Further- 
more, we can avoid the problem of constructing consistent, model-independent 
operators for conventional potentials by evolving to low momentum the po- 
tential and operators from chiral EFT. More generally, since Hartree-Fock is 
a reasonable starting point for many-body calculations, the large arsenal of 
techniques developed for the Coulomb many-body problem becomes available 
and should be explored as well. 
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